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Abstract 

Measurement of the three dimensional trajectory and specific ionization of recoil protons using a hydrogen gas time 
projection chamber provides directional information about incident fast neutrons. Here we demonstrate directional 
fast neutron detection using such a device. The wide field of view and excellent gamma rejection that are obtained 
suggest that this device is well suited to searches for special nuclear materials, among other applications. 
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1. Introduction 

Fast neutron imaging is a promising technique for 
Special Nuclear Material (SNM) search and safeguards 
applications due to the penetrating nature of fast neu- 
trons and since naturally occurring background rates 
are low |Q]]. Several techniques have been proposed 
for the directional detection of fast neutrons for these 
applications. These include Time Projection Cham- 
bers (TPCs) filled with H 2 [3), 3 He gas 0, or a 
mixture of gases J5). Considerable progress has also 
been made in adapting organic scintillator multiple pro- 
ton recoil imagers designed for solar neutron imag- 
ing |6| to the lower energies typical of fission neutron 
sources 0E|9). 

Here we describe the design, operation, and perfor- 
mance of a H2 filled TPC optimized for fast neutron de- 
tection. As in all neutron detectors [10], we observe the 
product of a nuclear reaction between an incident neu- 
tron and a nucleus in the target medium - in this case 
protons that have undergone elastic scattering. How- 
ever, not only can we detect the occurrence and energy 
of such a scatter, we can also track in three dimensions 
(3D) the direction and specific ionization of the recoil- 
ing proton. This novel capability provides significant in- 
formation about the incident fast neutron direction, es- 
pecially when averaged over several independent scat- 
tering events. 
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We demonstrate several features of this device that 
are of particular relevance to SNM search and safe- 
guards applications. First, the ability to localize sources 
at tens of meters standoff provides a clear benefit over 
non-imaging detectors, as well as providing an improve- 
ment in signal to background. Second, the TPC exhibits 
a wide field of view - in principle 4n, and in practice 
very nearly so - which is a feature of obvious benefit 
in a search application. Finally, the TPC possesses very 
good particle discrimination ability, particularly against 
gamma rays. 

2. Theory of Operation 

The TPC HUE) was invented in the 1970's and has 
been used extensively in fields as diverse as accelerator 
physics and dark matter searches. The TPC is capable 
of producing a full 3D reconstruction of the momentum 
and ionization produced by a charged particle in its ac- 
tive volume. Here we describe the application and opti- 
mization of this powerful detection technique to the task 
of fast neutron imaging. 

The charged particle tracking provided by TPCs can 
be used for directional neutron detection via several 
techniques. Tracking the charged reaction products of 
a capture reaction (e.g. 3 He(n,p) 3 H) yields the inci- 
dent neutron direction via the vector sum of the reaction 
product momenta - knowledge of the reaction 2 _va l ue 
also yields the incident energy. Similarly, recording the 
momentum vectors of two recoiling particles set in mo- 
tion via elastic scatters by the same neutron enables a 
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complete kinematic reconstruction of that neutron's in- 
cident direction and energy. Measurement of only one 
track and the vertex of the second recoil in a double elas- 
tic scatter constrains the incident neutron direction to 
lie upon a line. Finally, measurement of the momen- 
tum vector of a single recoiling particle from an elastic 
scatter provides a less complete constraint upon the in- 
cident neutron direction. In this situation, the tightest 
constraint is provided when the recoiling particle is a 
proton, because its mass is similar to that of a neutron: 
the incident neutron direction is constrained to lie in the 
hemisphere behind the recoil direction. Furthermore, 
it is worth noting that averaging an ensemble of single 
elastic proton scatters can fairly quickly determine the 
direction of a neutron source |3|. Given the relatively 
high efficiency for a single scatter occurring, this ap- 
proach appears attractive for rapid SNM discovery and 
search applications. This is the measurement mode that 
we have optimized for in this work. 

Among hydrogen bearing gases, alkanes possess the 
advantage of high hydrogen density as a function of 
gas pressure. This is advantageous in weight sensi- 
tive applications, since operation at a lower pressure 
allows a lighter pressure vessel, or indeed only a non- 
permable gas envelope at atmospheric pressure. Alka- 
nes are also commonly transported in pressurized form 
in commerce. 

H2 has previously been used in only a handful of 
TPCs fl31[l4l [B). Use of pure H 2 as the TPC gas for 
this application has the advantage that no carbon atoms 
are introduced into the TPC active volume. These in- 
crease the electron density, and therefore reduce track 
length making precise directional measurements more 
difficult, and provide large mass scattering centers that 
can result in large neutron deflections. For a fixed num- 
ber density of hydrogen atoms, H2 will result in the 
longest track and therefore provide the lowest energy 
threshold and most precise trajectory measurement. We 
have therefore chosen to use a gas primarily composed 
of H2 for this proof of principle demonstration. To en- 
sure good gas gain stability we use a gas mixture com- 
posed of 90% H 2 and 10% CH 4 (by volume). 

Of course, there are several other gases that could be 
used, but each has significant drawbacks. 4 He is readily 
available and non-flammable, but the average recoiling 
nucleus carries less energy and is not as well aligned 
with the neutron direction as is the case with a recoiling 
proton. Furthermore, the larger specific ionization of 
the recoiling 4 He nucleus would result in relatively short 
tracks. Finally, any radon or other alpha-decaying con- 
taminant in the counting gas or vessel materials could 
produce a high background indistinguishable from neu- 
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Figure 1: A schematic diagram of the TPC. 

tron induced recoils. 10 BF3 could be used in place of 
scarce 3 He, but it is toxic, and again the neutron capture 
reaction products are heavy, resulting in short, difficult 
to measure particle tracks. 

3. Detector Design 

A schematic diagram of the TPC system is shown 
in Fig. [T] A stainless steel gas vessel (50 cm length, 
46 cm diameter) encloses a field cage and charge read- 
out plane. To allow flexible operation of this prototype 
system, the gas vessel was designed for operation at 
pressures up to 10 atmospheres. To support this max- 
imum pressure the vessel wall is 6 mm thick, with end 
flanges of 3.8 cm thickness. As only non-mobile lab- 
oratory use was envisioned, no effort was expended to 
minimize the vessel mass. 

The field cage used to establish the drift field com- 
prises five FR-4 Printed Circuit Board (PCB) sheets 
with conductive strips etched upon them. One end is 
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maintained at a constant potential (-17.3 kV), while a 
resistive chain along the field cage establishes the drift 
field of 315 V/cm. The field cage defines an active gas 
volume of 25 cm x 25 cm x 40 cm = 25 liters. 

Ionization charge is drifted to an amplification and 
readout region that begins with a Frisch grid, main- 
tained at -3.6 kV. This is followed by a linear array of 
128 copper coated stainless steel wires of 30 /im diam- 
eter, orientated perpendicular to the horizontal plane of 
the TPC. A detailed study found that the TPC working 
gas would not embrittle these fine wires [2|. The wires 
are spaced at 2 mm intervals, and are held at ground po- 
tential. Behind the wire plane is PCB with 64 etched 
strips, each 4 mm apart, orientated parallel to the TPC 
horizontal plane. The strip plane is also maintained at 
-3.6 kV. 

The crossed wire and strip planes constitute the two 
dimensional charge readout of the TPC. Both wires and 
strips are connected to charge sensitive preamplifiers, 
directly in the case of the wires, and via de-coupling ca- 
pacitors for the strips. Due to a problem with the partic- 
ular connection scheme adopted during the data taking 
period reported here, only every other the strip could be 
read out, resulting in an effective 8 mm pitch. 

The TPC pressure vessel contains fitting for gas in- 
let and outlet, as well as pressure and temperature mea- 
surement. The system is operated in a constant pressure 
mode using a PID controller to regulate the position of 
a needle valve at the gas inlet. A port opposite the read- 
out plane allows for the introduction of a 242 Cm alpha 
particle source that is used for calibration studies. 

For the work described here, the TPC was operated 
at a pressure of 2 atmospheres, using a 90% H2/10% 
CH4 gas mixture. This operating pressure resulted 
in track lengths for 0.5 MeV recoil protons of about 
20 mm - well matched to the readout pitch and trig- 
gering scheme, and providing reasonable detection effi- 
ciency for an incident fission neutron spectrum (« 40% 
of recoils generated by an incident fission neutron spec- 
trum have an energy greater than 0.5 MeV). 

4. Data Acquisition 

Signals produced in the charge-sensitive preampli- 
fiers are digitized using a VME data acquisition system 
(DAQ) that is controlled using the ORCA data acquisi- 
tion software 07). The data acquisition hardware con- 
sists of nine 22-channel Waveform Digitizer/First Level 
Trigger (WFD/FLT) cards and a single Second Level 
Trigger/Interface card. These were originally produced 
by FZ Karlsruhe for the fluorescence telescopes of the 
AUGER experiment ifTBI . The WFD cards sample at 



10 MHz and record 1000 12 bit samples per event, i.e. 
record a total of 100 /-is per event. 

Each channel has an independent discriminator, 
which is applied to the average of 5 samples. The wave- 
forms of channels above threshold are recorded if 5 or 
more channels exceed threshold within 4 fis. The dis- 
criminator averaging, discriminator thresholds, trigger 
multiplicity and coincidence time can all be adjusted. 
The trigger parameter selections used here were found 
to be suitable for reducing spurious triggers on elec- 
tronic noise, while retaining a high efficiency for neu- 
tron interactions in the TPC. However, as a consequence 
of the coincidence requirement, short tracks at certain 
angles are not able to trigger data acquisition. This re- 
sults in a reduction in field of view of 10% (Sec. |5.5| l. 
We note that this is a consequence of the particular 
triggering scheme adopted for this proof of principle 
demonstration - a triggering and charge readout scheme 
sensitive to single channels, and the charge arrival time 
on that one channel would yield full 4n sensitivity. 

5. Event Reconstruction 

5.7. Waveform Conditioning 

Waveforms are first tested for ADC saturation, and 
the event rejected if it is found. The overall gain of 
the system is carefully adjusted so that neutron tracks 
do not cause such saturation - small discharges or a- 
particle tracks at steep angles with respect to the read- 
out plane are rejected by this selection. Next, a baseline 
subtraction is performed, using the average of the first 
and last 100 samples of each trace. An example of the 
waveforms that form a typical event are shown in Fig. [2] 

5.2. Hit Finding 

After baseline subtraction, a software threshold is ap- 
plied to each waveform to search for "hits" due to ion- 
ization tracks. Since the data set predominantly contains 
single elastic scatter events, we assume in this analysis 
that each waveform contains only one hit. The hit time 
is found by performing a quadratic fit over a 20 sam- 
ple range centered upon the maximum amplitude sam- 
ple in the waveform. A set of hits thus found for the 
wire and strip planes in a typical event are shown in 
Figures |3fa)&(b). 

The procedure described above finds the projection 
of the 3D ionization track onto the pair of crossed two 
dimensional (2D) readout planes. We are, of course, in- 
terested in reconstructing the full 3D track. The wave- 
form timing information is used to form coincidences 
between strip and wire plane hits. The waveforms from 
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Figure 2: The background subtracted waveforms of a typical event for 
(a) strips and (b) wires. 

each wire and strip containing a valid 2D hit are mul- 
tiplied together, and a hit search similar to that de- 
scribed above is performed, yielding a set of 3D hits 
(Fig. |3jc)). Waveforms that overlap in time produce 
combined traces with clear "hits"; those that do not 
overlap in time produce flat combined traces with no 
overlapping "hit". 

5.3. Tracking 

To determine track parameters, we begin with the as- 
sumption that proton recoil tracks are linear, i.e. will 
exhibit little straggling, which is appropriate given the 
light mass of the stopping medium. To identify linear 
features we apply a discrete Hough transform E2ll23l . 
since this technique is robust against noise. If a linear 
feature is found, hits that lie further than « 1 cm from it 
are discarded. 

We implement the Hough transform with fairly 
coarse bins for computational efficiency and operational 
simplicity. To realize the full tracking resolution of the 
TPC we then apply a three dimensional least squares 




Figure 3: An example event, comprising hits found in the (a) strip 
readout plane, the (b) wire readout plane, and (c) three dimensional 
hits deduced from those. 



fit to the remaining hit coordinates to extract the final 
(continuous) track parameters. These are expressed as 
two angles; an elevation angle referenced from the hor- 
izontal plane that runs through the center of the TPC 
(±90°, denoted "Elevation"), and an azimuth angle ref- 
erenced from the axis of the TPC vessel perpendicular 
to the charge readout plane (±180°, denoted "Phi"). Ex- 
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Figure 4: The ionization profile of the event shown in Fig. [5] 

amination of the track ionization profile determines the 
particle orientation along the track. 

5.4. Ionization Measurements 

With knowledge of the track parameters, the linear 
position of each hit along the track can be easily calcu- 
lated. The wire charge measurements are appropriately 
scaled to yield the specific ionization along the particle 
track (Fig. [4j. Examination of the track and identifica- 
tion of the position of the Bragg peak allows breaking 
of the two-fold pointing degeneracy. The point of great- 
est specific ionization is designated the Bragg peak, and 
the remaining ionization values are averaged to yield a 
quantity that we denote "Mean Tail Ionization." These 
two quantities are used to assess track quality and to 
perform particle identification. Integration of the spe- 
cific ionization along the length of the profile is used to 
determine the particle energy. 

5.5. Uniform Efficiency and Track Quality Cuts 

We examine several quantities to assess the quality of 
the tracks that are reconstructed, and to ensure a uniform 
angular acceptance. First, we exclude track angles and 
track lengths at which the data acquisition system is not 
uniformly efficient as a function of direction. In prac- 
tice, this means excluding track angles within 10° to the 
drift axis of the TPC (80 < |Phi| < 100), track angles 
at high elevations (|Elevation| > 80), and requiring that 
particle tracks be at least 20 mm long. The former re- 
quirement reduces the field of view by « 10%, while the 
later results in a loss of efficiency of about 25% when 
measuring fission neutron sources. We again note that 
this field of view reduction is a feature of this partic- 
ular prototype implementation, rather than of TPCs in 
general. 



Figure 5: The distribution of the Bragg Peak Location (BPL) parame- 
ter for a background run (red) and a run with a neutron source (black). 
Many background events do not have a Bragg peak at their terminus 
(BPL x or 1) and can therefore be rejected. 



We then examine the ionization profile of events to 
determine whether they follow the expected form of 
stopping charged particles, i.e. whether the event ion- 
ization profile resembles a Bragg curve. The location of 
the Bragg peak is determined, and the event rejected if 
it does not fall within the first or last 10% of the track 
length (Fig. [5j. 

5.6. Particle Identification 

Information from the ionization profile also allows 
identification of tracks created by protons, and rejection 
of gamma and alpha particle tracks (due to the inad- 
vertent use of a considerable area of lead solder on the 
interior of the field cage, there is a few Hz alpha par- 
ticle rate in TPC from 210 Pb decay). The first step of 
our particle identification technique is examination of 
Energy divided by Track Length (Average Specific Ion- 
ization). However, when short tracks dominated by the 
Bragg peak are considered, this quantity is a poor dis- 
criminant. This is demonstrated in Fig. [6ja), where data 
from neutron and alpha source calibration runs are com- 
bined, and compared to a prediction calculated using the 
SRIM package [ 19 1. At short length, the alpha and pro- 
ton bands converge. 

To provide robust particle discrimination at short 
track length, we also compare the Mean Tail Ioniza- 
tion quantity to the track length (Fig. |6|b)). For long 
tracks, this quantity is approximately independent of 
length. For shorter tracks, dominated by the rise in spe- 
cific ionization approaching the Bragg peak, the Mean 
Tail Ionization begins to rise, but distinct proton and al- 
pha bands remain. 
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Figure 6: The event track length is compared to (a) the particle energy 
and (b) Mean Tail Ionization, in order to distinguish particle type. In 
these plots data from separate neutron and alpha source calibrations 
runs are combined. The solid lines are predictions calculated using the 
SRIM package, while the dashed lines indicate the proton selection 
regions. 



A simple cut on the Mean Tail Ionization variable 
(40 < Mean Tail Ionization < 180) provides good dis- 
crimination (dashed lines in Fig.[6|b)), albeit at the ex- 
pense of some efficiency at short lengths. To provide a 
sense of average signal and background rates encoun- 
tered in this device, Fig.[7]compares the Mean Tail Ion- 
ization distribution of a background run and a run with 
a neutron source placed 8.7 m from the TPC. 



6. TPC Calibration 



A variety of measurements were made to characterize 
and calibrate the TPC in preparation for neutron mea- 
surements. 



Figure 7: The mean ionization found in the "tail" of tracks is used to 
identify particle type. Here we compare a run with a neutron source 
8.7 m from the TPC (black) to a background run of equal duration 
(red). The peak at x 100 is due to recoil protons while that at ss 250 
is due to alphas. 



6.1. Drift Speed Measurement 

Reconstruction of tracks in three dimensions re- 
quires measurement of the the speed at which ioniza- 
tion charge drifts to the readout plane. Our prototype 
TPC did not include a laser calibration system or exter- 
nal lifetime monitor; therefore the following calibration 
method was used. A neutron source was placed approx- 
imately 10 meters from the chamber at 10° increments 
in the horizontal plane. The drift speed loaded into the 
reconstruction software was then varied through reason- 
able values until the sources were measured to be at the 
correct locations. The drift speed depends on the pres- 
sure, temperature and gas composition but was typically 
around lcm/jus for the data presented here. The Mag- 
boltz program 1 18 1 was used to calculate the drift veloc- 
ity via numerical integration of the Boltzmann transport 
equation, yielding qualitative agreement with the exper- 
imental measurement. 

6.2. Channel-by- Channel Gain Correction 

The readout channel gains were calibrated relative to 
one another using neutron data. Event ionization pro- 
files, like that shown in Fig. [4] were fit to an analytic 
approximation of the form of the Bragg curve [21 1, 



dE 

dx 



(x) 



(1) 



where R is the range of the particle, a is a constant, and 
p — 1.8 for protons. The scale factor that would have 
to be applied to each point on the track to match the fit 
was calculated, and averaged over many events to cal- 
culate a gain correction factor. Typical corrections were 
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of magnitude 10 - 20%. This procedure was repeated 
iteratively several times, with the correction factor for 
one reference channel being held fixed. This calibra- 
tion was checked by applying a pulsed electrical signal 
directly to the Frisch grid. This yielded similar channel- 
by-channel corrections, implying that the major source 
of gain variation was variation among the preamplifiers. 

6.3. Ionization Loss due to Electronegative Contami- 
nants 

Outgassing from detectors components can introduce 
electronegative impurities that absorb primary ioniza- 
tion electrons as they drift to the charge readout plane. 
The effect of outgassing was measured by sealing the 
gas vessel and following the system gain over a pe- 
riod of 5 days. A gain degradation of almost 7%/day 
was observed. When a constant gas flow is maintained, 
an equilibrium is reached at which impurities are being 
added by outgassing at the same rate at which the gas 
flow is removing them. With a flow of 60 liters/hour 
it was found that this equilibrium was reached only 
12 hours after introducing the counting gas to an evac- 
uated TPC vessel, with a net reduction in gain of only 
4%. Therefore, at least 12 hours was allowed to elapse 
before any data taking. We note that no effort was made 
to select materials that exhibit low outgassing for this 
laboratory prototype. This suggests that use of materi- 
als that exhibit low outgassing, e.g. metals, kapton, and 
non-porous ceramics, could improve this situation con- 
siderably, and possibly allow for sealed operation. 

6.3.1. Correction for Gain Variation due to Tempera- 
ture Variation 
The small effect of ambient temperature variations on 
gain was also measured. Given that the TPC is operated 
in a constant pressure mode, temperature variations re- 
sult in gas density variations, which in turn result in gain 
variations 1201: 



(2) 



dG dp dT 
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Several resistive temperature devices mounted upon the 
strip plane were used to correlate changes in gas temper- 
ature with changes in gas gain over the course of several 
days. A linear relationship between temperature and 
gain was observed, with a coefficient of 1.5%/°C. Given 
that typical diurnal temperature variations in our labora- 
tory are as 0.2°C, this is a small effect, but nonetheless 
a temperature derived gain correction is applied. The 
application of thermal insulation to the exterior of the 
gas vessel would presumably reduce this small effect yet 
further. 
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Figure 8: The measured proton recoil spectrum due to a fission neu- 
tron source (black) compared to a Monte Carlo prediction (red). 



6.4. Energy calibration 

The conversion factor between the total measured 
ionization and energy for proton recoil events was de- 
termined using measured track lengths. The relation- 
ship between proton track length and proton energy was 
deduced for the gas composition and pressure used here 
using SRIM Ifl9l . The length/energy relationship thus 
determined was parametrized using a 2nd order polyno- 
mial, which was in turn fit to the proton data presented 
in Fig. |6ja) to yield the energy scale calibration in terms 
of MeV. The recoil proton spectra due to a 252 Cf fission 
neutron source can be seen in Fig. [8] This is compared 
to a simple Monte Carlo prediction that includes only 
those recoils that are fully contained within the active 
TPC volume. This condition effectively applies a high 
energy cutoff. 

7. TPC Design and Performance Simulation 

In order to assess the characteristics and pointing 
ability of our detector design, a Monte Carlo simula- 
tion of the TPC has been carried out using the MCNP 
(26) and MCNP-Polimi codes |27). The Monte Carlo 
model incorporates the main features of the detector de- 
sign (stainless steel chamber, field cage, and ionization 
gas), but not, in this implementation, laboratory walls 
or the intervening atmosphere. Simulation predictions 
of the angular resolution that can be achieved with the 
single elastic scatter technique are compared directly to 
data in Sec.lO 

We also used this simulation to predict the magnitude 
of scattering from the gas vessel and field cage. It is pre- 
dicted that ~ 92% of neutrons that scatter in the TPC ac- 
tive volume experience their first scatter in that volume, 
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with the remainder predominantly scattering first off the 
steel of the vessel. The small fraction of neutrons whose 
path is altered before interaction in the active volume is 
entirely acceptable in this proof of principle demonstra- 
tion, and could be reduced by careful tailoring of the 
gas vessel design to the exact parameters required by a 
particular application. 

8. Neutron Measurements 

8.1. Directional Response 

The imaging performance of the TPC was assessed 
by placing a 60 pCi 252 Cf fission neutron source at a 
variety of angular positions and standoffs (as large as 
17.7 m). In Fig. [9] we show the response of the TPC to 
three such experiments, demonstrating the broad field 
of view of this device. As discussed in 0, there is a di- 
rect correlation between the peak of angular distribution 
of recoil protons and the direction of incident neutrons. 
Thus, we bin the reconstructed direction of detected re- 
coil protons, and overlap a contour plot of the resulting 
distribution upon a panoramic photograph of the experi- 
mental laboratory. The correspondence of that distribu- 
tion with the source position can clearly be observed. 

To demonstrate the benefit of directional detection in 
increasing signal to background, consider the data pre- 



sented in Fig. 10 a), comparing data taken with the neu- 
tron source 17.7 m from the TPC and background. In 
this example, the excess of events due to the source 
is less than the background, when integrated over the 
entire angular range. However, the directional detec- 
tion technique demonstrated here yields a considerably 
greater signal/background ratio when that integral is 
limited to a narrower angular range centered upon the 
source direction. 

The angular resolution, i.e. ability to resolve closely 
spaced sources, of the simple source direction recon- 
struction method presented here (mean position of the 
recoil distribution) is just the width of that distribu- 
tion. We draw the readers attention to the apparent dif- 
ference in the width of the recoil distribution between 
the various experiments shown in Fig. [9] We tenta- 
tively attribute this to the effects of room scatter - in 
Figs.[9|b)&(c) the neutron source was close to concrete 
walls, cabinets, etc, while it was far from walls as was 
possible in Fig.|9ta). We therefore expect the data pre- 
sented in Fig. |9Fa) to be a better representation of the 
inherent angular resolution of the TPC, although further 
experimentation in a large, empty experimental location 
would be required to fully demonstrate this. This inter- 
pretation of the differences in recoil width is supported 
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Figure 10: (a) The recoil azimuth angle distribution of events acquired 
with a neutron source 17.7 m from the TPC (black), and background 
(red), (b) Background subtracted azimuth angle distribution of events 
due to neutrons from the 17.7 m distant source, compared to an MCNP 
prediction (red). 



by comparisons of the measured widths to that deter- 
mined from the MCNP simulation. This comparison is 



made in Fig. 10 b) for an experiment distant from walls, 



while that displayed in Fig. 1 1 was obtained with the 
neutron source adjacent a laboratory wall. 

While the recoil proton distribution itself is fairly 
broad, this is not an indication of the accuracy with 
which the mean direction of that distribution, i.e. the 
neutron source direction, can be determined. This is 



demonstrated in Fig. 12 where the direction within the 
TPC horizontal plane of a 252 Cf source is determined 
as a function of dwell time via a Gaussian fit to the re- 
coil distribution. As neutron counts are accumulated, 
the estimate of the mean position of the recoil distri- 
bution improves considerably. Detection of even a few 
tens of neutrons results in a significant reduction of the 
search space, compared to a non-imaging detector. 
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Figure 9: Contour plots demonstrating the pointing ability and wide angular field-of-view of the TPC with a 252 Cf source placed at various positions 
within the experimental laboratory. The blue star denotes the source position. The source standoffs were (a) 17.7 m, (b) 5.7 m, and (c) 4.0 m. Note 
that the TPC is senstive beyond Elevation = ±40; this representation of the data is limited by the abilities of the panoramic camera used. 



8.2. Neutron Detection Efficiency 

In Table [T] we examine event rates and the effect of 
various analysis steps in order to estimate the neutron 
detection efficiency of the TPC. Data collected with a 
252 Cf fission neutron source placed 8.7 m and 17.7 m 
from the TPC vessel are compared to background runs 
of equal length in each case. We use these data sets, 
since they are respectively likely to either overestimate 
(due to room scatter) or underestimate (due to several 
inches of unqualified intervening material) the effi- 
ciency. We estimate the number of neutrons produced 
by the source to be (3 + 0.3) x 10 s /s. The large un- 
certainty results from the advanced age of the source 
(24 years), and the consequent uncertainty as to the 



contribution of 25() Cf spontaneous fission to the neutron 
emission rate Il24l . Again, further experimentation in a 
large, empty experimental location will be required to 
gain a more accurate estimate. 

With no analysis cuts applied, an excess of events due 
to the neutron source equivalent to a a: 1 % neutron inter- 
action efficiency is observed. Application of all analysis 
cuts results in a » 50% reduction in efficiency. These 
experimental results are in reasonable agreement with 
calculations of the expected efficiency obtianed via both 
a first principle calculation and the MCNP simulation. 
These yield predictions for the raw neutron interaction 
efficiency of * 1.1% and for the ultimate efficiency after 
selection cuts of * 0.4%. 
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Table 1 : Estimates of the TPC neutron detection efficiency. Errors quoted reflect only counting statistics and neutron source strength uncertainty. 
The neutron flux reaching the TPC active area from the source at 8.7m distance is 32 ± 3 Hz and 8 ± 0.8 at 17.7 m, assuming no scattering or 
attenuation in the intervening material or atmosphere. 



Distance 


Analysis Cut 


Neutron + Bkgd Rate (Hz) 


Bkgd Rate (Hz) 


Neutron Rate (Hz) 


Efficiency (%) 


8.7 m 


Raw 


0.85 ± 0.02 


0.40 ± 0.02 


0.44 ± 0.03 


1.4 ±0.2 




Particle ID 


0.41 ± 0.02 


0.055 ± 0.006 


0.36 ± 0.02 


1.1 ±0.1 




Track Quality 


0.24 ± 0.01 


0.024 ± 0.004 


0.22 ± 0.01 


0.69 ± 0.08 




Uniform Efficiency 


0.19 ±0.01 


0.019 ± 0.003 


0.17 ±0.01 


0.55 ± 0.07 


17.7 m 


Raw 


0.377 ± 0.004 


0.350 ± 0.004 


0.027 ± 0.005 


0.35 ± 0.08 




Particle ID 


0.064 ± 0.002 


0.040 ± 0.001 


0.024 ± 0.002 


0.32 ± 0.04 




Track Quality 


0.033 ± 0.001 


0.019 ± 0.001 


0.015 ±0.001 


0.19 ±0.03 




Uniform Efficiency 


0.027 ± 0.001 


0.015 ± 0.001 


0.012 ±0.001 


0.16 ±0.02 



60 
50 
40 
30 
20 
10 




22 44 



Detected Neutrons 
176 340 



fi±+ + + r 



Mi 



-150 -100 -50 



50 100 150 
Phi (deg) 



Figure 11: Background subtracted azimuth angle distribution of 
events due to neutrons from a neutron source placed close to labo- 
ratory walls, compared to an MCNP prediction (red) 



We also note the low measured neutron background 
rate of 0.015 ± 0.001 Hz, which implies a background 
flux of ~ 8 x 10 -4 cm -2 s _1 , assuming the predicted de- 
tection efficiency of 0.4%. This is in reasonable agree- 
ment with quoted values for the integrated fast neutron 
background rate (« 1.5 x 10~ 3 cm~ 2 s _1 |Q~|) for fission 
neutron energy range (0.5 - 10 MeV), especially given 
the large uncertainties related to the neutron shielding 
provided by and cosmogenic production within the con- 
crete walls of the laboratory. 



8.3. Non-neutron Background Sensitivity 

To demonstrate the inherent gamma ray insensitivity 
of this device, we placed a 2 MBq 60 Co source beside 
the TPC. It was estimated that the interaction rate of 
60 Co gamma rays within the TPC gas was « 400 Hz. 
As can be seen in Fig. 13 no appreciable gamma re- 



sponse was measured. In the neutron selection region 
the event rate without the gamma source present was 




16 30 
Acquisition Time (min) 

Figure 12: The precision with which the direction of a neutron source 
located 8.7 m from the TPC can be determined improves as increased 
statistics are accumulated. The error bars correspond to the standard 
deviation of the azimuthal angle of the recoil protons for a source 
placed in a horizontal plane of the detector. 



0.022 ± 0.002 Hz, agreeing with that measured with the 
source present (0.020 ± 0.002 Hz). We therefore esti- 
mate the gamma rejection of this device to be > 10 6 . 

We have also examined the response of the TPC to 
microphonics, by striking the TPC vessel with a metal- 
lic object with enough force to set the readout wires in 
motion. Spurious events were generated, most likely 
by wire oscillation, but these did not pass the neutron 
selection cuts. The net result was a reduction is live- 
time, due to a saturation of the relatively low throughput 
DAQ system, e.g. tapping on vessel 6 times per minute 
resulted in a 25% livetime reduction. A higher through- 
put DAQ and a more rigid charge readout system, e.g. a 
LEM or thick GEM [25 1, would improve this consider- 
ably. Such effects are not observed during normal labo- 
ratory operation, e.g. personnel moving around the de- 
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Mean Tail Ionization (arb.) 



Figure 13: The mean ionization measured in the track tail is compared 
for a run with a 2 MBq 60 Co source beside the TPC (black) and a 
background run of equal duration (red). No appreciable gamma ray 
response is observed. The peak at ss 250 is due to alpha particle tracks 
within the TPC, while if they were measured compton electrons due 
to gamma rays would be expected to be less than 10. 

9. Conclusion 

In conclusion, a prototype directional detector for fast 
neutrons based on a H2 filled TPC has been designed, 
constructed, and characterized. The detector exhibits fa- 
vorable characteristics, such as the field of view nearly 
equal to the entire An solid angle, low sensitivity to 
gamma-ray background, and the ability to rapidly lo- 
calize neutron source directions. The ability of the TPC 
to point to a fission neutron source has been experimen- 
tally demonstrated in a laboratory environment with the 
efficiency of approximately 0.5%, and at several tens of 
meter standoff. 

Future efforts will focus on the exploration of higher 
resolution imaging capabilities using multiple neutron 
scatters and the incorporation of energy information into 
the reconstruction for single scatters, and on the opti- 
mization of detector design for field applications related 
to the detection of special nuclear material. In partic- 
ular, lightweight vessels operated at atmospheric pres- 
sure, using for example isobutane, seem attractive for 
mobile search applications. 
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